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A  NEW  METHOD  OF  ANALYSIS 
OF  COMPOSITE  STRUCTURES 


1  -  INTRODUCTION 


Composite  materials  consisting  of  high  tensile  resin-impre¬ 
gnated  fibers  are  being  more  and  more  frequently  used  in 
structures  capable  of  high  mechanical  performance.  Direct 
calculation  of  deformation  of  these  structures  using  the 
finite  elements  method  raises  major  difficulties  due  mainly 
to  the  very  high  number  of  heterogeneities  in  the  material. 
Computation  methods  are,  therefore,  based  on  investigation 
of  equivalent  homogeneous  materials,  i.e,  effective  behavior 
moduli  (Willis,  Hashin  ...). 

In  this  paper  we  use  the  homogenization  method.  This  me¬ 
thod  applies  when  the  material  being  investigated  has  a  pe¬ 
riodic  structure.  It  can  then  be  shown  that  when  the  dimen¬ 
sions  of  the  period  tend  homotheiically  to  zero  the  fields  of 
deformation  and  stresses  tend  to  those  corresponding  to  a 
homogeneous  structure  whose  elastic  properties  can  be 
computed  precisely  when  a  single  period  of  the  composite 
medium  to  be  investigated  Is  known.  This  boundary  value 
structure  is  the  homogenized  struaure  and  its  behavior 
coefficients  are  the  homogenized  coef '!':ients.  This  is  the 
macroscopic  equivalent  structure.  Furthermore  by  a  locali¬ 
zation  procedure  the  method  allows  an  easy  computation 
of  the  microscopic  field  of  stresses  and,  in  particular,  of 
stress  forces  at  the  boundaries  between  fibers  and  matrix. 
These  stress-forces  are  particularly  important  because  they 
can  initiate  cracks  and  delarninations.  The  overstresses  at 
the  microscopic  level  may  produce  fiber  ruptures. 


)  Material  reinforced  by  periodically  arranged,  parallel 
fibers  (Figure  1) 


f^ig,  1  :  PARALLEL  FIBERS 

ii)  Material  consisting  of  a  very  large  number  of  parallel 
layers  of  homogeneous  materials  superposed  periodi¬ 
cally  (Figure  2) 


After  presenting  the  general  method  of  homogenization,  Fig,  2  :  MULTIPLE  LAYERS 
which  leads  first  to  an  homogenized  equivalent  macroscopic 
structure  and  secondly  to  a  localization  procedure  for  com¬ 
puting  the  field  of  microscopic  stresses  and  stress  forces,  v/e  This  is  followed  by  the  numerical  results  obtained  by  using 

apply  the  method  to  two  types  of  composite  materials  :  the  MODULEF  code. 
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2-  DESCRIPTION  OF  THE  HOMOGENIZATION 
METHOD  (1)  [4]  [51  llOl  l12l 


2.1  —  Formulation  of  the  problem 


Let  US  consider  an  elastic  body  which  occupies  a  region  Q 
related  to  a  system  of  orthonormal  axesOx-j  X2  x^.  This 
body  is  subjected  to  a  system  of  voluminal  forces]  fj  |  and 
surface  forces  |  Fj  [  on  a  portion  fj:  of  boundary o  The 
other  portion  of  the  boundary  is  to  which  a  zero  move¬ 
ment  condition  is  imposed. 


The  field  of  stresses  at  equilibrium  satisfies  the  equilibrium 
equations 

1 )  +  f.  =  0  in  n 

2)  (Tjj  n.  =  F,  on  r 

Furthermore,  the  material  is  elastic  with  fine  period’C  struc¬ 
ture,  i.e.  n  is  covered  by  a  set  of  identical  periods  of  rec¬ 
tangular  (Fig.  4)  or  hexagonal  (Fig.  5)  or  more  complicated 
shape  such  as  the  examples 


given  in  Figures  6  and  7. 


Fig.  6 :  MATERIALS  WITH  FINE  PERIODIC  STRUCTURE 


Fig.  7 :  MA  TE RIALS  WITH  FINE  PERIODIC  STRUCTURE 

All  the  period  forms  must  be  such  that  opposing  faces  which 
correspond  in  a  translation  can  be  defined  two  by  two. 

In  all  cases  we  shall  designate  as  Y  a  period  characteristic 
of  the  material  which  has  been  enlarged  by  homothetics 
and  fixed  once  and  for  all.  e  then  designates  the  homothetic 
ratio  which  is  small  and  which  takes  us  from  Y  to  a  period 
in  the  elastic  material.  The  elastic  structure  of  the  material 
is  then  fully  i^nown  if  it  is  given  over  a  single  period,  e.g.the 
enlarged  period  Y  related  to  the  onhonormal  axis  system 
Oy^  ^2  ^3*  ®ijkh  ^7^  the  coefficients  of  elasticity 

on  Y,  which  generally  alter  very  quickly  with  respect  to  y, 
but  satisfy  in  all  respects  the  symmetry  relation 

®ijkh  ^  ®jikh  ~  ®khij 
and  positivity  relation 

3^0  >  (y)  Tjj  Tj.  T. 


Fig.  4  :  MATERIALS  WITH  FINE  PERIODIC  STRUCTURE 


Fig.  5  :  MATERIALS  WITH  FINE  PERIODIC  STRUCTURE 


The  functions  y  — (y)  defined  on  Y  are  extended  by 
Y’periodicity  to  the  entire  space  Oy^  ^2  Vs  assumed  to  be 
covered  by  contiguous  periods  identical  to  Y. 

The  coefficients  of  elasticity  in  the  material  H  are  then 
^ijkh  ^  defined  by  : 

®iikh  =  =ijkh<Y>-  Y  =-f 


For  greater  simplification  in  the  text  we  shall  write  : 

3  (vl  =  {  Sijkh  (y>  i  .  a  f  (xl  =  a  (-f-  I  .  0  =  I  0;j  I 

and  shall  consider  a  ly)  or  a  ^  (x)  as  known  matrix  6x6 
indexed  by  the  symmetrical  pairs  (i,  j). 


The  law  of  elasticity 

j 

3>  %•  =  ‘Sikh  *  “kh 


is  written 

as  a'(x)  e(u). 


where  .  t,,  £ 


When  an  ambiguity  is  possible,  either  (u)  or  (u)  will 
be  specified  depending  on  whether  the  drift  occurs  with 
respect  to  x  or  y.  The  boundary  conditions  are  finalized 
by 

4)  u  =  0  on 

The  problem  posed  by  ID  (2)  (3)  (4)  has  a  unique  solution 
which  depends  on  6  and  which  we  shall  designate  u  ®  ;  to 
this  corresponds  a  field  of  stresses  O  ^  given  by  : 

5)  r  (x)  e  (u 

Numerically  it  is  very  difficult  when  e  is  small  to  calculate 
u  ^  since  there  are  a  large  number  of  heterogeneities  in  the 
elastic  medium.  We  therefore  try  to  obtain  a  limited  expan¬ 
sion  of  the  solution 
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0  ®  {x,  y)  =  s  (y)  By  (u®) 

or '  (x,  y)  =  8  (y)  (u'l  +  (u®lj 

o  ^  (x,  y)  s  8  (y)  ^By  (u^)  +  8^  (uh] 

The  equilibrium  equations  (1)  applied  to  <7^  give 
a  *  1-  +  f  -  =  0 

or  in  a  more  condensed  form 

9)  div  +  f  =  0  . 

Given  the  expansion  (8)  of  we  have  (^*) 

10)  -V  div^  +  div  a^) 

^  e  ^ 

+  div  y  -f  div  x  ^  *"  “ 

xC53ay6:Y. 

The  boundary  conditions  {2)  are  treated  in  the  same  way 

ID  —  a®  .  n  +  .n  -  F  +  €  0^  .n  +  =  0 

€ 

for  X  e  Tp  y  e  Y. 


2.2  Asymptotic  expansions 

The  solution  is  affected  by  two  factors  : 

i)  The  first  is  the  scale  of  H  and  arises  from  the  forces 
applied  and  the  conditions  at  the  boundaries. 

ii)  The  second  is  due  to  the  periodic  structure  ;  it  is  on 
the  same  scale  as  the  period  and  is  repeated  periodical¬ 
ly. 


Finally  the  conditions  (4)  mean  that 

12)  u  ^  +  e  u^  +  ...  =  0 

forxerQ,yeY. 

By  making  the  various  powers  of  €  zero  we  obtain  : 


div  y  a  ®  =  0 
a  0  =  a  (V)  ey  (u  ®) 


This  justifies  looking  for  an  asymptotic  expansion  of  the 

'  ( 

[  d'''  V/ 

0  ^  + 

div 

o 

II 

o 

o 

: 

14)^ 

1 

y 

X 

form  : 

la’  = 

aly) 

[®y 

(u’)  +  e^(u°)] 

6)  u  ^  =  u  ®  (x,  y)  +  6  (x,  y)  +  6 ^  u^  (x,  y)  4-... 

div 

1  V 

0^  4- 

div 

^  a  ’  +  f  =  0 

where  the  u^rlx,  y)  are,  for  each  xGH  ,  Y  periodic  func¬ 

15)| 

1  y 

a  ly) 

[*V 

(u  +  e  lu  ' )1 

tions  with  respea  to  the  variable  y  G  Y.  Then  y  =  is  ap- 

L  y 

*  J 

plied  to  (6).  Associated  with  the  expansion  (6)  is  an  expan¬ 
sion  of  the  field  of  deformation  e  (u  (*). 

7)  e  (u  ey  lu^l  +  (u^)  +  e^  (u^ ) 

+  e  lu^)  -f-  By  (u^)]  4-  ... 

and  of  the  field  of  stresses  o  ^ 

8)  O  ^  =  -1  ^  y)  +  ^  ^  y^  e  0  ^  lx,  y)  4-  .. 


The  equations  (11)  and  112)  will  be  used  later. 

(*)  Note  that 

-i-  u“(x,Y):-f^  u*U,yl  +  1 u®(x,y) 

clx|  0X|  e  e.y| 


t  - 
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2  J  -  Reiolutioni  23J  div  j,  <  >  +  f  =  0  in 


The  $v«em$  (13)  (14)  (15)  contain  differential  operators 
in  y.  They  therefore  constitute  equations  wit)'  partial 
derivatives  on  the  period  of  base  Y,  the  unknown  factors 
being  the  Y-periodic  functions. 

System  (13)  :  This  leads  immpdiately  to  : 

16)  a®=0,u®=u®  (x) 


If  we  introduce  Z  z  <a  .we  have 


24)^ 


[div^  I  +f  :0 

^  •!  ^  jj 


’kh 


(u 


Using  equation  (12)  and  taking 
obtain  : 


in  n 


the  mean  on  Y  in  (11),  we 


System  (14)  :  In  view  of  (’6)  it  is  reduced  to  : 

17)  diVy  a’ r  0,  o’ =  a(y)  [ey  (u ’)  +  e^(u°)j 

The  deformation  {\P)  is  a  function  only  of  x  ;  it  there¬ 
fore  plays  the  role  of  a  parameter  with  respect  to  the  diffe¬ 
rential  system  in  y.  Due  to  the  linearity.  CJ*  ,  may  there¬ 
fore  be  written  in  the  form  : 


18) 


^1  _  jkh  Q  Oj 


M  ^  (y)  (u 


where 


'  kh 


(U°)=±(£A  +  ^ 


2  5 


div =  0 


19)<s’'^=  a(y)[c'^^  +  By  {X^^\] 
(  is  Y'periodic 
The  tensor  has  components  given  by 


=  0 

\z.nzf 


on  , 
on  Tp  . 


The  system  (24)  with  boundary  conditions  (25)  is  a  well 
posed  elasticity  problem  ;  the  equilibrium  equations  are 
unchanged,  as  well  as  the  boundary  conditions.  The  elastic 
constitutive  relation  is 


-  kh 


kh 


(u^ 


It  is  homogeneous  since  the  coefficients  q  ij  given  by 
(22)  are  independent  of  .  These  coefficients  define 

the  equivalent  homogeneous  material.  They  are  called  homo¬ 
genized  coefficients.  The  stress  field  Z  s:(-S  jj)  is  called  the 
macroscopic  stress  field  and  is  defined  by 

I  =  <a’> 

T)ie  strain  field  E  z  6,^  (u  °)  is  called  the  macroscopic 
Strain  field  and  satisfies 

E  z  <e^  (u°)  +  By  (u’)> 


c''*’  =7<‘^ik^jh  +  ^ih  ^jk' 

ij  •<: 

It  can  be  proved  that  the  system  (19)  determines  the  vec¬ 
tor  X^^  (yl  to  within  an  additive  constant. 

For  any  function  4)=  <f>  (x,  y  I,  we  define 

Jy 

The  solution  a  ^  of  (14)  is  given  by, 


Ich 

It  can  be  proved  that  the  homogenized  coefficients 
satisfy 

f  ='»kh  ‘='’iikh> 

(3ai>o,  sjj  s,^h^a,  Sjj  $ij.  VsjjZSji 

This  shows  that  (q  )  are  reasonable  elastic  coefficients 
and  that  the  macroscopic  scale  problem  (24)  (25)  has  a 
unioue  solution. 


20)  0  ’  (x,y)  z  a  (y)  -  ey  a''^)]  e  kh 

and  taking  the  mean  value,  we  obtain, 

1  kh  ,  Os 

21)  <o.'i>z  q  e  Ku  (u 

'I  i| 

where 

22)  =<^ikh'v)>-<aiipq(v)  ep^  (A''''’(y)> 

System  (15)  :  It  suffices  to  take  the  mean  on  Y  in  the 
the  first  equation  to  obtain 


23  —  Microscopic  fields.  Localization 

The  stress  field  o  ^  (x,  y)  is  the  first  term  of  the  asymptotic 
expansion  (3)  of  the  stress  field  a  ^  (x)  solution  of  the  ini¬ 
tial  exact  problem.  The  field  o  ^  {x,  y|  is  called  the  micros¬ 
copic  stress  field.  If  we  imagine  that  at  each  point  xd2, 
there  is  a  small  FY  period  with  Its  composite  structure,  then 
0  ^  (x.  yl  gives,  for  x  kept  fixed  in  H  .  a  stress  field  in  this 
period. 


88-4 


It  can  be  shown  that  o  ^  (x)  -  o  Mx  ,  )  tends  to  zero 

in  the  L  '  (fi)  norm  when  e  tends  to  zero.  This  proves 
that  0  Mx.  I  is  a  good  approximation  of  a  ^  (x)  when  6 
is  small.  The  microscopic  stress  field  a  ^  (x,  y)  ,  y  r  —can 
be  calculated  as  follows :  ^ 

i)  First  we  obtain  the  six  ..V  (y)  vector  fields  on  Y, 

each  one  been  associated  with  tensor  These 

six  vector*fietds  are  solution  of  problem  (19)/which  is 
an  elastic  type  problem  on  the  inhomogeneous  period 
Y. 

ii)  From  the  vector  fields  (y)  we  get  the  homogeni¬ 
zed  coefficients  q  by  formula  (22). 

>1 
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Calculation  of  the  homogenized  coefficients  calls  for 
the  resolution  of  (19).  In  the  present  case  the  coefficients 
aiikh  (y)  are  independent  of  73  ;  the  result  is  that  the  fields 
A ‘I  (y)  are  also  independent  of  y3  ;  in  (19)  the  varipus 
indices  give  a  zero  contribution  when  they  refer  to 
making  computation  of  A*'i  (y)  a  bidimensional  problem.  ^ 

3.2  —  Numerical  results 

In  ail  the  cases  studied,  the  homogenized  material  is  ortho- 
tropic,  in  other  words  the  law  of  behavior  has  numerous 
zero  elements  as  shown  in  the  table  below  : 


iii)  We  solve  the  macroscopic  scale,  homogenized  elastic 
problem  (24)  (25)  on  Q  .  \t  gives  the  macroscopic 
stress  field  Z  (x)  and  the  maaoscopic  strain  field 
ej^(u^)  z  E(x),  for  xSn. 

iv)  Localization  procedure  :  using  formula  (20)  we  can 

calculate  o  ‘  (x,  y).  For  x  fixed  in  ,  this  stress 
field  on  Y  shows  how  the  macroscopic  stress 
Z(x)=:  (x,  y)>  is  localized  in  an  e  Y  period  at 

xeQ 

It  can  be  proved  that  when  e  tends  to  zero,  the  stress  field 
a  ^  (x)  tends  to  X  (x)  in  the  weak  L^  iQ)  topology.  Never¬ 
theless  a  ^  (x,-~)  is  a  better  approximation  of  a  ^  (x)  than 
Z  (x)  :  the  norm  L  ^  iCl)  convergence  implies  that 

(7  ^  (x)  -  a  ^  (x,^) 

tends  to  zero  for  almost  every  point  in  H  ,  while  the  weak 
(Hi  convergence  does  not.  The  macroscopic  stress  field 
Z  (xl  is  just  a  mean  value  while  takes  into  ac¬ 

count  the  fine  periodic  structure  of  the  composite  material. 

3  -  APPLICATION  TO  AN  ELASTIC  MATERIAL  REIN¬ 
FORCED  BY  FIBERS  RUNNING  IN  THE  SAME 
DIRECTION  (11(7) 


3.1  —  Principle 


\ 

^1111  ^1122  ^1133  ®  0  0  j 

j 

i 

^2222  **2233  0  0  0 

1^22/ 

(  - 

^3333  0  0  0 

X 

r33\ 

( 

SYM.  2c2323  0  0 

] 

r23( 

] 

2q,3.3  0 

^^13  j 

/ 

2q,2^2 

1 

1 

^^12/ 

where  { cr  -  ^  and  |  ^  jj  [  are  stress  and  strain  tensors. 

The  law  (26)  is  inverted  conventionally  to  be  written  [6]  : 


The  computations  of  the  previous  paragraph  are  applied  to 
an  elastic  material  formed  from  a  multitude  of  resin-impre¬ 
gnated  unidirectional  fibers  whose  geometric  distribution  is 
periodic  in  a  plane  perpendicular  to  their  direction  X3, 


(») 
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bringing  out  the  following  : 

The  Young's  moduli  E.j.  E2.E3  in  the  directions  of  ortho- 
tropY 

The  Poisson's  coefficient  ^23' *^13' ^12 
The  shear  moduli  G23.  0^3.  G.j2 

The  numerical  results  which  follow  have  been  obtained  by 
using  the  MOD'JLEF  code  (2)  .  They  have  been  produced 
for  numerous  values  of  the  ratio  of  impregnation  and  various 
forms  of  the  cross  section  of  the  fibers. 

We  give  her©  a  part  of  the  results  obtained  for  various  forms 
of  fiber,  and  also  the  curves  showing  the  change  in  these 
coefficients  with  respect  to  the  ratio  of  resin  impregnation 
for  fibers  of  circular  section  (Fig.  9-10-11). 
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aligned  circular  fiber 

FIBER 

El  =  3.8  10®  MPa  623  =  2.10^  MPa 
E2=E3z:.145  10®  G^2=Ql3  =  3.8  lO'^MPa 

^12=^13  =.22  U23='25 

resin 

E::  3520  MPa  i;  =  ,38 
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Fig,  W:  VARIATION  OF  TRANSVERSE  YOUNG 
MODULI 
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Fig.  9  :  VARIA  T/ON  OF  LONGITUDINAL  YOUNG's 
MODULI 


Fig,  1 1  :  VARIA  TION  OF  SHEAR  MODULI 
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-  Anisotropy  curves  (Fig.  11) 

It  is  impoaant  to  note  that  the  homogenized  media  obtained 
are  generally  not  transversally  isotropic.  This  comment  Is 
clearly  demonstrated  if  the  Young's  modulus  is  calculated  in  a 
transverse  direction  with  polar  angle  d.  By  applying  the 
Young's  modulus  on  vector  radius  we  obtain  the  curves 
given  in  Figure  12.  For  the  material  to  be  transversally 
isotropic,  the  curves  plotted  should  be  arcs  of  a  circle  cen¬ 
tered  at  the  origin. 

Note  : 

The  Young's  modulus  in  direction  9  is  given  by  : 


+  sin  ^9  cos  ^  5  ^ 


h 


^23 


This  relation  enabled  the  anisotropy  curves  in  Figure  12  to 
be  plotted. 

The  material  is  transversally  isotropic  if  E  (0)  is  not  depen¬ 
dent  on  0,  which  is  equivalent  to 


E,r  2023(1  -fK, 


RESIN  IMPREGNATION  RATIO  BY  VOLUME  50  % 
(FIBERS //TO  XI). 


FIBER 

E1=  380000  MPa  ;  G23  =  20000  MPa  ;>'23::.25 
E2z  14500  MPa  ;G13- 38000  MPa  ;>*13=.22 
E3rE2  ;G12  =  G13  ;Kl2i>'l3 


RESIN 
E=3520  MPa 

Ki:.38 


E> 

"  1  « 

ri7 

,n  1  C23 

012  1  C13 

192009 

9730 

9070 

23 

2« 

.30 

2452 

5597 

3334 

191500 

8290 

•  100 

40 

29 

.30 

2623 

*315 

3347 

0 

191500 

1 

7620 

! 

7620 

,44 

29 

29 

2755 

3062 

3662 

Fig.  12  :  TRANSVERSE  ANISOTROPY  FOR  2  CROSS 
SECTIONS  OF  FIBER 

3.4  —  Stagger 

If  the  fibers  are  staggered,  i.e.  if  a  period  characterizing  the 
material  has  the  form  shown  in  Figure  13,  we  obtain  diverse 
characteristics  in  accordance  with  the  relative  values  of  the 
sides  of  lengths  of  the  rectangular  cell. 

i)  If  1  r  1  (square  cell)  :  the  characteristics  of  directions 
0y2  and  Oy^  are  identical,  and  have  the  same  Young's 
modulus  in  particular. 

ii)  If  Ji  z  V3,  i.e.  if  the  fibers  are  located  at  the  apexes  of 
an  equilateral  triangle  (Fig.  13)  it  can  be  shown  that 
tile  material  is  transversally  isotropic.  This  property  is 
true  for  any  impregnation  level  of  the  resin. 

iii)  The  bisecting  direaions  0  y2  and  0  ^3  play  the  same 
roles  irrespective  of  the  values  of  land  the  impregna-  • 
tion.  In  particular,  the  Young's  moduli  and  £3  in  • 
these  directions  are  always  equal. 

iv|  In  Figure  14  are  plotted  the  Young's  and  shear  moduli 
corresponding  to  the  various  values  of  Jl  varying  from 
1  to  2  and  for  the  same  resin  impregnation  level  by 
volume.  ForJ,=1,  the  cell  is  square  and  naturally  = 

^2  E  ^  £2  foil-iYS  since  then  the  fibers 

are  at  the  apexes  of  an  equilateral  triangle  and  the  ma¬ 
terial  is  then  transversaily  isotropic,  which  implies 
In  the  same  figure  are  plotted  the  values 


!i  = 


J2- 


E^  z  ^2*^^  Young's  modulus  in  the  bisector  direc¬ 
tions  0  y  ^  and  0  72- 

For  1  z  \/5  we  find  a  triple  point  since  naturally  the 
transverse  isotropy  then  implies 


E,=  E2=E^  =  E2 


i 
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CHARACTERISTICS 
FIBER  :  E  84000  MPA 
.22 

RESIN  :  E  4000  MPA 
.34 

RESIN  RATIO  36 
FIBERS //Y1 
REF.  1  :0Yi  Y2  Y3 
REF.2:0Yi  Y2  Y3 


3.5  —  Comparison  with  experiments 

The  development  of  this  previsional  method  is  aimed  at  ob¬ 
taining  complete  sets  of  characteristics  for  three-dimensional 
computations  of  composite  structures  through  the  finite 
elements  method. 

The  possibilities  of  experimental  characterization  are  indeed 
very  reduced.  Few  tests  are  reliable,  each  one  being  specific 
to  a  characteristic,  not  permitting  to  reach  them  all.  The 
results  of  measurements  being  very  scattered  in  relation  to 
production  batches,  mean  values  have  to  be  used. 


The  extreme  variety  of  resins  give  a  very  wide  range  of  pro¬ 
ducts  to  be  used  in  production.  Each  fiber-resin  pair  can  be 
associated  within  variable  proportions.  It  is  unthinkable  to 
be  able  to  experiment  all  configurations. 

Each  material  is  therefore  charaaerized  in  an  incomplete, 
dissimilar  and  inaccurate  manner. 

Tables  presented  hereafter  explain  application  of  the  homo¬ 
genization  theory  to  the  two  materials  :  glass  R  -  Resin 
Ciba  920  (36  %  -  Resin  in  volume)  and  carbon  CTS  -  Resin 
Qba  920  (50  %  resin  in  volume).  We  have  considered  several 
distributions  and  shapes  of  fiber. 


Fig,  13  :  EQUIDISTANT  STAGGER  (L::1.73)  Taking  these  values  into  account,  average  measured  values 

were  assigned  to  glass-resin  composites  while  values  obtained 


Fig.  14  :  VARIATION  OF  YOUNG'S  AND  SHEAR  MODULI 
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by  transposition  of  tests  results  and  proportion  computa¬ 
tion  were  assigned  to  carbon-resin  composite.  Asa  reminder, 
characteristics  obtained  with  two  bidimensional  previsional 
methods  :  PUCK  [11]  and  HALPIN-TSAI  (14]  were  also 
given.  For  reasons  indicated  formerly,  comparisons  must  be 
cautiously  made.  Results  obtained  for  glass-resin  composite 
with  staggered  fibers  layout  at  the  apexes  of  an  equilateral 
triangle  (ensuring  transverse  isotropy)  are  nearest  to  measu¬ 
red  values  .With  those  two  methods,  )^23  ^23 

obtained. 

As  far  as  carbon  based  composite  is  concerned,  it  is  less 
clear  but,  in  this  case,  the  real  shape  of  the  fiber  is  not  ob¬ 
served.  On  the  other  hand,  when  the  shape  is  more  accurate 
(t Kidneys  shaped),  the  direaion  of  the  fiber  does  not  vary 
and  is  therefore  as  little  realistic.  Of  course,  a  configuration 
taking  into  consideration  random  direction  will  probably  be 
nearer  to  the  truth. 

For  the  two  considered  materials,  estimates  based  on  the 
homogenization  theory  are  nearer  to  those  based  on  the 
widely  used  HALPIN-TSAI  method. 

The  homogenization  theory  seems  efficient  to  compute  the 
mechanical  characteristics  of  composite  materials. 


Validity  of  the  results  is  evidently  subjected  to  the  assump¬ 
tions  made  on  shapes  and  lay-out  of  fibers.  However,  the 
undeniable  advantage  of  this  method  aims  at  supplying 
complete  and  consistent  sets  of  values,  mutually  coherent. 


COMPARATIVE  TABLE  FOR  CARBON  CTS 
RESIN  CIBA  920  (50  %  RESIN  IN  VOLUME)  COMPOSITE 


REFERENCE 

VALUES 

H0MOGENIZAT10f4  THEORY 

OTHER  PREVISIONAL  METHODS 

ALIGNED 

CIRCULAR 

FIBERS 

STAGGERED 

CIRCULAR 

FIBERS 

•KIDNEY" 

SHAPED  FIBERS 

(MEAN  VALUES) 

PUCK 

- 

HALPIN-TSAI 

El 

(  MPa) 

120  000 

119  299 

119  293 

, 

119  290 

113  260 

119  250 

E2 

(MPa) 

6  000 

6  284 

li  035 

8  000 

11  620 

5  620 

E3  {  6  OOQ 

(MPa)  1 

6  284 

6  035 

7  950 

11  620 

5  623 

— 

Yi2  i  0.28 

0.299 

1 

0.299 

0.31 

0.3 

0.3 

Yi3 

0.28 

0.299 

0.299 

0.29 

1 

0.3  1  0.3 

.  1 

Y23  j  0.20 

0.435 

0.457 

0.27 

1 

-  !  - 

Gi2 

(  MPa ) 

3  800 

3  454 

3  391 

4  500 

i 

4  250  {  3  350 

Gi3 

(  MPa) 

3o00 

3  454 

3  391 

3  200 

4  250  j  3  350 

/  MP  *  1 

■Ql 

3  266 

2  100 

1 

-  i 

(MPa) 


3  2G6 


2  100 


•  • 


ORiG^K":.  y-i 

COMPARATIVE  TABLE  FOR  GLASS  R-RESIN  OP  OiJALriY 

'«ia  A  non AT  DCCiM  tiLj  %/ni  I  iiiACi  r^rtiiO/^ciTC  V 


CIBA  920  (36  %  RESIN  IN  VOLUME)  COMPOSITE 


MEASURED 

VALUES 


HOMOGENIZATION  THEORY 


align:  : 


OTHER  PREVISIONAL  METHODS 


El 

(MPal 

55  000 

E2 

(MPa) 

17  000 

E3 

(MPa) 

17  000 

>'12 

0.26 

y^2 

0.26 

>'23 

- 

Gi2 

(MPd) 

5  600 

Gi3 

(MPa) 

5  600 

G23 

(MPa) 

— 

20  275 
(^  =  13  498) 


STAGGERED 
CIRCULAR  FIBERS 

PUCK 

HALPIN-TSAI 

55  215 

54  450 

54  450 

16  016  . 

18  800 

18  570 

16  016 

18  800 

18  570 

0.256 

0.264 

0.264 

0.229 

(>^3=0.487) 

a357 

6  383 

5  337 

6  383 

5  387 

4  539 

(G'23z8  250) 

5  832 

3.6  —  Microscopic  stn*«  field 

Given  a  structure  consisting  of  a  unidirectional  material  and 
subjected  to  a  simple  shearing  overall  stress  field  within  the 
plane  (1.2)  normal  to  the  direction  of  fibers,  the  biaxial 
stress  tensor  at  macroscopic  level  is  : 


0  o  0 
a  0  0 
0  0  0 


The  localization  method  allows  calculation  of  both  the  stress 
field  at  microscopic  level  which,  in  any  point  of  the  tratertal 
period,  is  : 


Fig,  15  :  STRESS  FORCES 


^12  ^ 
a-f2  022  ^ 

0  0  On 


and  the  stress  forces  at  fibre  /  matrix  interface  as  represented 
in  Figuka  15. 
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4  -  Aff  LICATK)N  TO  A  PERIODIC  STACK  OF 
HOMOOEWeOUS  LAYERS  [4]  [7j 

4.1  -  Principi* 

We  shall  consider  a  periodic  stack  of  a  multitude  of  homo* 
genized  layers.  Each  layer  is  characterized  by  a  direction  of 
the  fibers.  In  the  stack  these  directions  vary  periodically 
whilst  remaining  orthogonal  to  axis  Ox^. 

interface: 

^iikh  '*3>  /  laytr  (p)=‘='’ijkh  (“"swnt) 


Fig.  16 !  MUL  TtPLF  LA  YFRS.  EACH  LA  YER  POSSESSES 
A  PLANe  OF  ELASTIC  SYMMETR  Y  NORMA  L 
TO  THE j(j AXIS  (MONOCLINIC SYMMETRY) 


In  this  situation  the  homogenization  formulae  are  conside¬ 
rably  simplified  since  the  problem  (19)  is  then  reduced  to  a 
system  of  differential  equations  which  may  be  solved  ex¬ 
plicitly.  For  the  details,  refer  to  D.  Begis,  G.  Duvaut,  A. 
Hassim  (1 1  arvi  to  the  references  in  this  publication. 


HOaiOGENIZED  MODULI 
OF  EACH  LAYER 

HOMOGEKlZl-O  MODULI 
OF  COMPOSITE 

El 

120  QC3  MPa 

45  128  MPa 

E2 

S  OCO  MPa 

45  128  MPa 

E3 

6  000  MPa 

6  198  MPa 

K23 

0.20 

0.188 

>•13 

0.28 

1  0.188 

1  _  . 

>•12 

0.28 

1  0.30 

G23  * 

2  500  MPa 

3  015  MPa 

G13 

3  800  MP. 

3  015  MPa 

G12 

3  800  MPa 

17  290  MPa 

Conclusion 

We  have  presented  several  applications  of  the  homogeniza¬ 
tion  techniques  for  computing  the  coefficients  of  elasticity 
of  composite  materials.  Other  applications  using  the  locali¬ 
zation  procedure  are  contemplated  as  regards  fine  analysis 
of  the  field  of  stresses  using  asymptotic  expansions,  the 
effect  of  defects  in  the  composites  (9l  and  more  generally, 
damage  to  the  materials  of  composite  structure  containing 
inclusions  or  precipitates. 

Strictly  speaking,  these  techniques  apply  only  to  absolutely 
periodic  structures,  but  with  the  backing  of  statistical  ana¬ 
lyses  it  is  possible  to  identify  the  fluctuations  likely  to  be 
produced  by  periodic  defects.  It  is  noted  generally  thai 
strict  periodicity  reinforces  the  anisotropy  of  the  computed 
homogenized  material  with  respect  to  the  industrial  model. 


4J2  —  Numeric*!  application 
As  an  illustratior^  vw  consider  two  cases  : 

1)  a  laminate  cor-s-sting  of  3  identical  layers  disposed  pe¬ 
riodically,  The  layers  have  equal  thickness  and  their 
fibers  orientat  ons  are  respectively  ~  60^.  0°,  and  60^ 
with  respect  to  the  -  axis.  The  homogenized  material 
then  prevsms  a  transverse  isotropy  which  complies 
with  the  ^e.'-ieral  results  on  isotropy,  cf.  (8)  . 

2)  a  laminate  c:ors.sting  of  18  layers  identical  to  the  above 
and  laid  yp  at  successive  angles  of  10^  to  each  other.  It 
is  to  checked  that  the  same  result  is  obtained  as  in 
the  previous  case. 

We  give  in  the  tafcie  presented  hereafter  the  moduli  of  each 
layer  and  the  moduli  of  the  composite  which  are  identical  6 
for  the  two  (3  layers  and  18  layers). 
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